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We present a variational many-body wave function for repelling bosons in rotating traps, focusing 
on rotational frequencies that do not lead to restriction to the lowest Landau level. This wave 
function incorporates correlations beyond the Gross-Pitaevskii (GP) mean field approximation, and 
it describes rotating boson molecules (RBMs) made of localized bosons that form polygonal-ring-like 
crystalline patterns in their intrinsic frame of reference. The RBMs exhibit characteristic periodic 
dependencies of the ground-state angular momenta on the number of bosons in the polygonal rings. 
For small numbers of neutral bosons, the REM ground-state energies are found to be always lower 
than those of the corresponding GP solutions, in particular in the regime of GP vortex formation. 
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Recent experimental advances in the field of trapped 
ultracold neutral bosonic gases have enabled control of 
the strength of interatomic interactions over wide ranges 
P, 13, 0, J from the very weak to the very strong. This 
control is essential for experimental realizations of novel 
states of matter beyond the well known Bose-Einstein 
condensate |3H,0. In this context, the linear ID Tonks- 
Girardeau regime of impenetrable trapped bosons has 
generated intensive theoretical activity |3, |^ and several 
experimental realizations of it have been reported most 
recently [ll^. 

Here we address the properties of strongly-repelling 
impenetrable bosons in rotating ring-shaped or 2D har- 
monic traps. To this end, we recall that impenetra- 
ble bosons are "localized" relative to each other 0, Q 
and exhibit nontrivial intrinsic crystalline correlations 
0. For a small number of bosons, N, these crystalline 
arrangements are reminiscent of the structures exibited 
by the well-studied rotating electron molecules (REMs) 
in quantum dots under high magnetic fields Conse- 
quently, we use in the following the term rotating boson 
molecules (RBMs). A central result of our study is that 
the point-group symmetries of the intrinsic crystalline 
structures give rise to characteristic regular patterns (see 
below) in the ground-state spectra and associated angu- 
lar momenta of the RBMs as a function of the rotational 
frequency for neutral bosons (or the magnetic field for 
charged bosons). 

An unexpected result of our studies is that the rotation 
of repelling bosons (even those interacting weakly) does 
not necessarily lead to formation of vortices, as is famil- 
iar from the case of rotating Bose-Einstein condenstates 
(BECs). In particular, for small N, we will show that the 
Gross-Pitaevskii energies (including those corresponding 
to formation of vortices) remain always higher compared 
to the ground-state energies of the RBMs. Of course, we 
expect that the rotating BEG will become the preferred 
ground state for sufficiently large N in the case of weakly 
repelling neutral bosons. We anticipate, however, that it 
will be feasible to test our unexpected results for small N 



by using rotating optical lattices, where it is established 
that a small finite number of atoms can be trapped per 
given site 01 • 

In a non-rotating trap, it is natural to describe a lo- 
calized boson (at a position Rj) by a simple displaced 
gaussian 0. When the rotation of the trap is consid- 
ered, the gaussian needs to be modified by a phase factor, 
determined through the analogy between the one-boson 
Hamiltonian in the rotating frame of reference and the 
planal motion of a charged particle under the influence of 
a perpendicular magnetic field B (described in the sym- 
metric gauge). That is, the single-particle wave function 
of a localized boson is 
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with Q = z/(2A^) and the width of the Gaussian A is 
a variational parameter; A = Is = \/hc/{eB) for the 
case of a perpendicular magnetic field B, and A = = 
^/Ti/ (2mri) in the case of a rotating trap with rotational 
frequency $1. Note that we consider a 2D trap, so that 
r = {x,y) and R = {X,Y). The hamiltonian corre- 
sponding to the single-particle kinetic energy is given by 
(r) = (P — ?i,Q X r)^/ {2m), for the case of a magnetic 
field, and by Hk{y) = (P - ?iQ x vf/{2m) - 
for the case of a rotating frame of reference 

A toroidal trap with radius tq can be specified by the 
confining potential 



(2) 



with Iq — y/h/ (mwo) being the characteristic length of 
the 2D trap. For n ^ 2 and Iq/tq this potential 
approaches the limit of a toroidal trap with zero width, 
which has been considered often in previous theoretical 
studies (see, e.g., ReL ^^). In the following, we consider 
the case with n = 2, which is more realistic from the 
experimental point of view. In this case, in the limit 
ro = 0, one recovers a harmonic trapping potential. 
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projected ground-state energy is given by 
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FIG. 1: Properties of A'^ = 8 neutral repelling bosons in a 
rotating toroidal trap as a function of the reduced rotational 
frequency Q/uio- The confining potential is given by Eq. Q 
with n = 2 and radius ro = 3lo, and the interaction-strength 
parameter was chosen as Rs = 50. (a) RBM ground-state 
energies, E^^'' . The inset shows the range < Q,/ljo < 0.3. 
The numbers denote ground-state magic angular momenta, 
(b) Energy difference E^^-' - _bUbhf ^-^^ r^^^^^j ^j^g^iaj. 

menta associated with (i) the RBM ground states [thick solid 
line (showing steps and marked as PR J); online black] and 
(ii) the UBHF solutions (thin solid line; online red) 



To construct an RBM variational many-body wave 
function describing N impenetrable bosons in the 
toroidal trap, we use N displaced orbitals ip{r,'R.i), i = 
1, 2, A'' [see Eq. (Q] centered at the vertices of a reg- 
ular polygon. Then, we first construct an unrestricted 
Bose Hartree-Fock (UBHF) permanent l^"™""") oc 
'/'i(rii)¥'2(ri2)---'/'iv(i"jiv)- The UBHF permanent 
breaks the circular symmetry of the many-body hamil- 
tonian. The "symmetry dilemma" is resolved through 
a subsequent "symmetry-restoration" step accomplished 
via projection techniques |^ 0, i-^-: we construct 
a many-body wave function with good total angular 
momentum by applying the projection operator Vl = 
(l/27r) /g de exp[i6'(L - L)], so that the final RBM wave 
function is given by 
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is the original UBHF permanent rotated by 
an azimuthal angle 9. We note that, in addition to hav- 
ing good angular momenta, the projected (PRJ) wave 
function j^*^^) has also a lower energy than that of 
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' (0)) and n{e) = 
the latter term ensures 
proper normalization. 

The many-body hamiltonian is given by 7Y = 
j:f=i[HKiri) + y(r,)] + E^, Vir,,r,), with the inter- 
particle interaction being given by a contact potential 
Vs — g6{Ti — Tj) for neutral bosons and a Coulomb po- 
tential Vc = Z'^e^/\Ti — Yj\ for charged bosons. The 
parameter that controls the strength of the interparti- 
cle repulsion relative to the zero-point kinetic energy is 
given by Rs ~ gmj (2tt11^) B for a contact potential and 
Rw = Z^e^/{^olo) foi" a Coulomb repulsion. 

For a given value of the dimensionless rotational fre- 
quency, fl/oJo, the projection yields wave functions and 
energies for a whole rotational band comprising many an- 
gular momenta. In the following, we focus on the ground- 
state wave function (and corresponding angular momen- 
tum and energy) associated with the lowest energy in the 
band. 

Fig. 1(a) displays the ground-state energy -Eprj of 
N = 8 bosons in a toroidal trap as a function of the 
dimensionless rotational frequency fl/oJo, with uiq being 
the trap frequency. The prominent features in Fig. 1(a) 
are: (i) the energy diminishes as Q/ujo increases; this is 
an effect of the centrifugal force, and (ii) the i?pRj curve 
consists of linear segments, each one associated with a 
given angular momentum L. Most remarkable is the reg- 
ular variation of the values of L with a constant step of 
N units (here A'^ = 8) [see inset in Fig. 1(a) and Fig. 
1(c)]. These preferred angular momenta L = kN with 
integer k, are reminiscent of the so called "magic angu- 
lar momenta" familiar from studies of electrons under 
high-magnetic fields in 2D semiconductor quantum dots 
[a 0. The preferred angular momenta reflect the in- 
trinsic molecular structure of the localized impenetrable 
bosons. We note, that the (0,8) polygonal-ring arrange- 
ment is obvious in the single-particle density associated 
with the UBHF permanent [see Fig. 2(b)]; (0,8) denotes 
no particles in the inner ring and 8 particles in the outer 
one. After restoration of symmetry, however, the single- 
particle (sp) density is circularly symmetric [see the PRJ 
sp-density in Fig. 2(c)] and the intrinsic crystallinity be- 
comes "hidden"; it can, however, be revealed via the 
conditional probability distribution 0,11 [CPD, see Fig. 
2(d)]. We note the Gross-Pitaevskii sp-density in Fig. 
2(a), which is clearly different from the PRJ density in 
Fig. 2(c). 

The internal structure for charged bosons in a toroidal 
trap (not shown) is similar to that of neutral bosons (Fig. 
2), i.e., a (0,8) ring arrangement, portrayed also in the 
stepwise variation (in steps of 8 units) of the total angular 
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FIG. 2: Single-particle densities and CPDs for N — 8 bosons 
in a rotating toroidal trap with fl/uo = 0.2 and Rs = 50. The 
remaining trap parameters are as in Fig. 1. (a) GP sp-density. 
(b) UBHF sp-density exhibiting breaking of the circular sym- 
metry, (c) RBM sp-density exhibiting circular symmetry, (d) 
CPD for the RBM wave function [PRJ wave function, see Eq. 
©I revealing the hidden point-group symmetry in the intrin- 
sic frame of reference. The observation point is denoted by 
a white dot. The RBM ground-state angular momentum is 
Lz — 16. Lengths in units of Iq. The vertical scale is the same 
for (b), (c), and (d), but different for (a). 



momenta. The internal structure is also reflected in the 
variation of the ground-state total energy as a function 
of the magnetic field. In contrast to the case of neu- 
tral bosons, however, the ground-state energy curve for 
charged bosons is not composed of linear segments, but of 
intersecting inverted-parabola-type pieces; this is due to 
the positive contribution of the Lorentz force compared 
to the negative contribution of the centrifugal force in a 
rotating trap. 

For RBMs in rotating harmonic traps, the polygonal- 
ring pattern of localized bosons becomes more complex 
than the simple (0, A'^) arrangement that appears natu- 
rally in a toroidal trap. Indeed, in harmonic traps, one 
anticipates the emergence of concentric ring structures. 
For = 6 neutral bosons in a harmonic trap, we ob- 
serve that, as in the case of a toroidal trap, the ground- 
state energy as a function of the reduced rotational fre- 
quency, ri/wo, [Fig. 3(a)] is composed of linear segments, 
but now the corresponding magic angular momenta [Fig. 
3(b)] vary in steps of — 1 = 5 units. This indicates 
an RBM consisting of two polygonal rings; denoted as a 
(1, 5) structure, with the inner ring having a single boson 
and the outer ring five. 

In Fig. 4(a), we display the RBM and mean-field GP 
ground-state energies of TV = 6 strongly repelling (i.e., 
Rs = 50) neutral bosons in a harmonic trap as a func- 
tion of the reduced angular frequency of the trap. The 
GP curve (thin solid line; online red) remains well above 
the RBM curve (thick solid line; online green) in the 
whole range < U/uiq < 1. The RBM ground-state 
angular momenta exhibit again the periodicity in steps 
of five units [Fig. 4(b)]. As expected, the GP total an- 
gular momenta are quantized [Lz = (no-vortex) or 
Lz = 6 (one central vortex)] only for an initial range 
<n/ujQ < 0.42. For Q/loq > 0.42, the GP total angu- 



lar momentum takes non-integer values and ceases to be a 
good quantum number, reflecting the broken-symmetry 
character of the associated mean field, with each kink 
signaling the appearance of a different vortex pattern of 
p-fold symmetry (p = 1,2, 3, 4,...) 0]; see an example in 
Fig. 4(c). 

The energetic superiority of the RBM wave function 
over the GP solution demonstrated in Fig. 4(a) was to 
be expected, since we considered the case of strongly re- 
pelling bosons. Unexpectedly, however, for a small num- 
ber of neutral bosons the energetic advantage of the RBM 
persists even for weakly repelling bosons, as illustrated 
in Fig. 5(a). Indeed, Fig. 5(a) displays the RBM (thick 
solid line; online green) and GP (thin solid line; online 
red) ground-state energies for A'^ = 6 neutral bosons in a 
trap rotating with ^I/luq = 0.85 as a function of the inter- 
action parameter Rg. The surprising result in Fig. 5(a) 
is that the GP energy remains above the RBM curve 
even for Rg 0. Of course the RBM wave function 
is very close to that of a BEC without vortices when 
Rs ^ (BECs without vortices are approximately fea- 
sible for small A^). However, for small N, our results 
show that BECs with vortices (i.e., for Lz > N) are not 
the preferred many-body ground states; instead, forma- 
tion of RBMs is favored. Note that the energy difference 
E'^^ — E^^'^ increases rapidly with increasing Rs, reflect- 
ing the fact that the RBM energies saturate (as is to be 
expected from general arguments), while the GP ener- 
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FIG. 3: Properties of A'^ = 6 neutral bosons in a rotating 
harmonic trap as a function of the reduced rotational fre- 
quency Q/ujo- The confining potential is given by Eq. Q 
with n — 2 and ro — 0, and the interaction-strength parame- 
ter was chosen as Kg = 50. The intrinsic molecular structure 
is (1,5). (a) RBM ground-state energies, E^^'. The inset 
shows a smaller range. The numbers denote ground-state an- 
gular momenta, (b) Total angular momenta associated with 
(i) the RBM ground states (thick solid line showing steps; on- 
line black) and (ii) the UBHF solutions (thin solid line; online 
red). 
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FIG. 4: Properties of GP solutions (thin solid line; online 
red) versus those of RBM wave functions (thick solid line; 
online green) for A'' = 6 neutral bosons as a function of the 
reduced rotational frequency fl/uuo- A harmonic trap is con- 
sidered, and the interaction strength equals Rs — 50. (a) 
Ground-state energies, (b) Associated ground-state angular 
momenta, (c) GP (BEG) sp-density at Q/uJo — 0.65 having 7 
vortices with a 6-fold symmetry (thus exhibiting breaking of 
the circular symmetry), (d) RBM sp-density at fl/cuo ~ 0.65 
which does not break the circular symmetry, (e) GPD of the 
RBM at Q/ujo = 0.65 revealing the intrinsic (1,5) crystalline 
pattern. The white dot denotes the observation point tq. 
Note the dramatic difference in spatial extent between the 
GP and RBM wave functions [compare (c) with (d) and (e). 
Lengths in units of lo- The vertical scale is the same for (d) 
and (e), but different for (c). 



gies (even with vortices fully accounted for) exhibit an 
unphysical divergence as i?^ ^ oo [Fig. 5(a)]; we have 
checked this trend up to values of Rs = 100 (not shown). 
Of interest again is the different behavior of the RBM 
and GP ground state angular momenta [Fig. 5(b)] (see 
also discussion of Fig. 4(b)). 

In conclusion, we have introduced (and studied the 
ground-state properties of) a variational many-body 
wave function for repelling bosons in rotating traps that 
incorporates correlations beyond the Gross-Pitaevskii 
mean-field approximations. This variational wave func- 
tion describes rotating boson molecules, i.e., localized 
bosons arranged in polygonal-ring-type patterns in their 
intrinsic frame of reference. For small numbers of neutral 
bosons, and in particular in the case of GP vortex for- 
mation, the RBM ground-state energies are lower than 
those associated with the corresponding Gross-Pitaevskii 
BEC solutions. Given the large differences between the 
properties of the RBM and BEC wave functions (which 
become more pronounced for larger interaction parame- 
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FIG. 5: Properties of GP solutions (thin solid line; online 
red) versus those of RBM wave functions (thick solid line; 
online green) for N = 6 bosons as a function of the interaction 
strength Rs- A harmonic trap is considered, and the reduced 
rotational frequency equals Q/ujo — 0.85. (a) Ground-state 
energies (b) Associated ground-state angular momenta. 



ter Rs), and the recently demonstrated ability to exper- 
imentally control Rs Q, IM Q ; we anticipate that our 
results could be tested in experiments involving rotat- 
ing optical lattices. Detection of RBMs could be based 
on a variety of approaches 0], such as the measure- 
ment of the spatial extent [contrast the RBM and BEC 
spatial extents in Figs. 4(c)-4(e )], o r the use of Hanbury 
Brown-Twiss-type experiments |l7l| to directly detect the 
intrinsic crystalline structure of the RBM. 
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